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The ground-state fidelity has been introduced recently as a tool to investigate quantum phase 
transitions. Here, we apply this concept in the context of a crossover problem. Specifically, we 
calculate the fidelity susceptibility for the BCS ground-state wave function, when the intensity of 
the fermionic attraction is varied from weak to strong in an interacting Fermi system, through the 
BCS-BEC crossover. Results are presented for contact and finite-range attractive potentials and for 
both continuum and lattice models. We conclude that the fidelity susceptibility can be useful also 
in the context of crossover problems. 
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I. INTRODUCTION 

A quantum phase transition is an abrupt change of the 
ground state of a many-body system when a controlling 
parameter A of the Hamiltonian crosses a critical value 
A c . It is then quite natural to expect that the overlap 
F(X + 5X, A) = |(#(A + <5A)|*(A))| between the ground 
states corresponding to two slightly different values of 
the parameter A, should manifest an abrupt drop when 
the small variation SX crosses A c . Such overlap, which 
has been named "ground-state fidelity" in the literature, 
should thus provide a tracer of a quantum phase transi- 
tion. An appealing feature of the ground-state fidelity is 
that it does not rely on the explicit knowledge, or even 
the very existence, of an order parameter associated with 
the quantum phase transition. 

The ground-state fidelity approach to quantum 
phase-transitions has been tested in a variety of 
modeb 1 ! 2 ! 3 ! 4 !^ 6 ^^^ 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 , and has been 
found to be quite effective in signalling the presence of 
a quantum phase transition. In particular, even phase 
transitions which are not associated directly with a lo- 
cal order parameter, like transitions topological in nature 
or of the Beresinskii-Kosterlitz-Thouless type have been 
found to be detectable by the ground-state fidelity (or 
derived quantities like the fidelity susceptibility) 3 - ^ 13 ' 19 . 
[Some phase transitions of high order in terms of the con- 
trolling parameter A have been found, however, to escape 
the analysis based on the ground-state fidelity-^.] 

In the present paper we make a step further in the 
study of the potentialities of the ground-state fidelity, 
by analyzing its applicability beyond the domain for 
which it was originally conceived. We consider specifi- 
cally a crossover problem, that is, a many-body system 
for which a substantial change in the nature of its ground 
state occurs over a finite range of the controlling param- 
eter A, rather than abruptly at a critical value A c . On 
physical grounds, we expect that the sudden drop of the 
ground-state fidelity at a critical value A c should in this 
case be replaced by a minimum, located in the region 
of A where the ground state is changing more rapidly, 
namely, in the crossover region. In terms of the fidelity 
susceptibility (defined shortly below), the divergence at 



a critical point should correspondingly be replaced by a 
peak, whose width should be associated with the width 
of the crossover region. 

To test the above ideas we consider specifically the 
BCS-BEC crossover, namely, the evolution of the ground 
state of a fermionic system in the presence of an at- 
tractive potential which is progressively increased in its 
strength. This evolution, first studied in Refs. I20ll2lll22l . 
has been studied then quite extensively, firstly in con- 
nection with the physics of high-temperature cuprate 
superconductor a 23 ' 24 ' 25 ' 26 ' 27 ' 28 , and more recently with 
ultracold Fermi atoms, where this crossover has been re- 
alized experimentall y 29 ' 30 ' 31 ' 32 ' 33 ' 34 . Several models sup- 
porting this crossover will be analyzed, both in con- 
tinuum and discrete space. Such a systematic analysis 
will allow us to draw conclusions of sufficient generality, 
which we believe to be applicable also to other crossover 
problems. 

The paper is organized as follows. In Sec. II we in- 
troduce the fidelity susceptibility and derive its expres- 
sion for the BCS wave- function. This quantity is calcu- 
lated explicitly in Sec. Ill for the BCS-BEC crossover in 
continuuum models. We consider specifically the three- 
dimensional contact and Nozieres-Schmitt-Rink poten- 
tials, and compare the information extracted from the 
fidelity susceptibilty with what is already known in the 
literature for the BCS-BEC crossover in these models. 
A similar analysis is presented in Sec. IV for the attrac- 
tive Hubbard model (in both two and three dimensions) . 
Section V gives our conclusions. 



II. THE FIDELITY SUSCEPTIBILITY FOR THE 
BCS WAVE FUNCTION 



The ground-state fidelity F(X + SX, A) depends, by its 
definition, on both the controlling parameter A and its 
variation, 5X. The somewhat artificial dependence on 
the actual value of the parameter 8X can be eliminated by 
considering the limiting expression for the ground-state 
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fidelity when 8X approaches zero. For small S A: 
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where the state |\&(A)) is assumed to be real and normal- 
ized. A sudden drop of the ground-state fidelity at the 
critical point will then correspond to a divergence of the 
fidelity susceptibility 2 -^: 
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Note that, in order to deal with a meaningful quantity in 
the thermodynamic limit, the expressions on the right- 
hand side of Eqs. ([3]) and ([J| in the definition of the 
fidelity susceptibility have been divided by the system 
volume f2. For a sufficiently large system, one has in 
fact lnF(A, A') oc f2. Indeed, barring the case when a 
correlation lenght £ diverges (i.e. when A or A' sit ex- 
actly at a critical point), the system can be thought as 
a collection (tensor product) of many identical subsys- 
tems of volume L D (with L 3> £ and where D is the 
spatial dimension). The overlap between two different 
ground states |^(A)) and |^(A')) will be then the prod- 
uct of the overlaps in each individual subsytem, yielding 
F(X,X') cx f(X, X') N ", where /(A, A') is the overlap in 
each subsystem and N s = fl/(L D ) is the total number 
of subsystems. This implies then lnF(A, A') oc f2 except 
when A or A' are exactly at a critical points. 

In this paper, we are interested in calculating the fi- 
delity susceptibility x(A) across the BCS-BEC crossover. 
Previous studie s 21 ' 22 ' 35 have shown that the BCS wave- 
function provides a reasonably good approximation for 
the ground-state wave-function over the whole BCS-BEC 
crossover, from the weak-coupling limit of highly overlap- 
ping Cooper pairs to the strong-coupling limit of tightly 
bound dilute composite bosons. 

We will thus calculate x(A) for the BCS wave-function: 
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Here, creates a fermion in the single-particle state of 
wave- vector k, spin a and energy ek, Uk and Vk are the 
usual BCS coherence factor v 2 = 1 — u\ = (1 — ^/ Ej.)/2, 



with £ fc = e fc - (j,, E k = + A 2 ,, where \x is the chemi- 
cal potential and A& the BCS gap function. The param- 
eter A in Eq.© stands generically for the appropriate 
coupling strength of the attractive interaction V\(k, k') 
which is driving the BCS-BEC crossover. The depen- 
dence of Uk and Vk on A in Eq.© is determined by the 



gap function and chemical potential fi, which depend 
on A through the gap and particle number equations: 



A* = - 



dk ' W) Afe ' 



(2tt) d 
dk 



2Ek> 



{2-kY 



■2vi 



(6) 
(7) 



where n is the particle number density. 

When the BCS wave function is inserted in Eq. ([J| for 
X(A) one obtains: 
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which, after some manipulations, can be written 



X(A) 



dk 



dfi dA k 
Ak dX + ^ k ^X 



(8) 



(9) 



In the next two sections we will analyze the behaviour 
of the fidelity susceptibility for several type of attrac- 
tive interaction V\{k, k'), by solving the coupled Eqs. © 
and ([7]) and by calculating then xW as determined by 
Eq. ^j. Section III will consider continuum models, 
while section IV will deal with a lattice model (the at- 
tractive Hubbard model). In that case the integration 
over k in Eqs. ©,© and © will be limited to the first 
Brillouin zone. 



III. CONTINUUM MODELS: CONTACT AND 
FINITE-RANGE POTENTIALS 

The simplest model Hamiltonian for the BCS-BEC 
crossover describes a system of fermions of mass m in 
continuum space, mutually interacting via a contact (5- 
like) interaction: 
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+ g /dr^(r)V{(r)V|(r)VT(r). (10) 



The Hamiltonian (fT0|) leads to ultraviolet divergencies, 
as it can be seen in the gap equation ^ when both the 
gap function and the interaction do not depend on wave 
vector (as it occurs for a contact potential) . These diver- 
gencies are, however, eliminated by expressing physical 
quantities in terms of the two-body scattering length a p 
rather than the bare coupling g. 

The above Hamiltonian (with the appropriate ultravi- 
olet regularization) has been studied quite extensively in 
the context of the BCS-BEC crossover, especially after 
the advent of experiments with ultracold Fermi atoms in 
the presence of a Fano-Feshbach resonanc e) 29 ' 30 ' 31 ' 32 ' 33 ' 34 . 
For these systems the Hamiltonian (ITUl) can in fact be de- 
rived from first principles 36 , as the effective Hamiltonian 
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FIG. 1: Fidelity susceptibility (in units of kp 3 ) for the three- 
dimensional contact potential as a function of the dimension- 
less coupling strength (SiFttF) -1 - 



describing the physics of the relevant degrees of freedom 
close to the Fano-Feshbach resonance. 

For a three dimensional contact potential the gap equa- 
tion ([6]) , when expressed in terms of the scattering length 
ap reads: 

m fdk.fl m\ 
~ 4^f~ = J J2nf \2E^ ~ WJ ■ 

The coupled Eqs. and determine the gap A and 
chemical potential /i in terms of the scattering length 
Of or, better, of the dimcnsionless coupling parame- 
ter (fcfaf) -1 , which is normally used as the coupling 
strength parameter for the 3D contact potential (here 
kp = (S^n) 1 / 3 , such that the scattering length ap is 
compared with the average interparticle spacing kp 1 ). 

In terms of this parameter the BCS and BEC limits 
correspond in principle to the conditions (/cf^f) 1 <C 
— 1 and (kpap)^ 1 ^> 1, respectively. Previous 
studie o 37 ' 38 i 39 ' 40 i 41 have shown, however, that the 
crossover between the above two different physical sit- 
uations is limited in practice to the rather narrow region 
-1 < (kpap)- 1 < 1. 

We have calculated the fidelity susceptibility © , where 
as a controlling parameter A we have taken the di- 
mensionless coupling strength (i.e. A = (kpap)~ r ), 
by solving the coupled Eqs. (fTTjl and Q (which can 
be suitably expressed in terms of elliptic integrals^) 
to determine A(A), fi(X) and then x(A) through the 
BCS-BEC crossover. The resulting fidelity susceptibil- 
ity x((kpap) _1 ) is presented in Fig. 1. 

The fidelity susceptibility presents a rather symmet- 
ric peak, which is located in the middle of the crossover 
region (specifically, the peak position corresponds to 
(kpap)^ 1 ~ —0.14) and whose half-height width marks 
precisely the borders of the crossover region —1 < 
(kpap)^ 1 < 1 mentioned above. Note that the size of 



the crossover region for the three-dimensional contact 
potential was determined in previous studies by calculat- 
ing specific physical quantitities, like e.g. the chemical 
potential, the BCS gap, or the superfluid critical tem- 
perature and by comparing their numerical values with 
analytic expressions valid in the BCS and BEC limits, 
respectively. This empirical way of defining the range of 
the crossover region, even though physically sound, could 
be criticized because of some degree of arbitrariness in 
choosing the physical quantity to look after, or the quan- 
titative criterion to conclude that a specific asymptotic 
(BEC or BCS) expression has been effectively reached. 

The plot of the fidelity susceptibility ^(/cf^f) -1 ) of 
Fig. 1 provides a somewhat more "intrinsic" criterion to 
locate the position and width of the crossover region, 
since it is not based on a specific physical quantity but 
on the measure of the rapidity of change of the ground- 
state wave-function through the crossover. The agree- 
ment between the position and width of the crossover 
region, as determined by the fidelity susceptibility, with 
previous results obtained with more empirical criteria to 
characterize the crossover region corroborates these pre- 
vious results while proving, at the same time, the utility 
of the ground-state fidelity for studying also crossover 
problems. 

We pass now to consider a finite-range potential. 
Specifically, we consider the separable potential intro- 
duced by Nozieres and Schmitt-Rink^ (NSR): 



y/i + (k/k Q ) yi + WW 

where k$ sets the range of the potential in momen- 
tum space. The finite range of the NSR potential al- 
lows for the occurrence of the density induced BCS-BEC 
crossover— which is instead not possible in the case of a 
contact potential. 

We have calculated the fidelity susceptibility for the 
NSR potential (with A = V) for various values of the 
particle density (which can be parametrized by the ratio 
between A:f and the momentum range of the potential 
fco). For each density, our calculated x(V) are peaked at 
a value V = V p whose position indicates where the rate 
of change of the BCS wave-function with respect to V is 
maximal. As already argued above, the value V p should 
thus be located in the crossover region. 

Fig. 2 compares the peak position V p extracted from 
the fidelity susceptibility for several values of kp/ko 
(full line) with the curves defined by the conditions 
fc-F^pair = 27T (dashed line) and fcpCpair = l/ir (dash- 
dotted line), obtained previously in Ref. 42) and intro- 
duced there to characterize the width of the crossover 
region (thus defining a sort of "phase diagram" for the 
BCS-BEC crossover). In particular, £ pa j r represents the 
pair correlation length, as defined in Ref. l26l and pro- 
vides an estimate of the Cooper pair radius. The BCS 
region is characterized by large overlapping Cooper pairs 
(such that fcF^pair 3> 1) while for the BEC region, with 
small nonoverlapping boson, /cF^pair 1- The two values 
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FIG. 2: Crossover "phase diagram" for the three-dimensional 
NSR potential. The full curve was obtained by determining 
the position of the peak in the fidelity susceptibility x(V) a t 
given values of the parameter Jcf /feo ■ The dashed and dash- 
dotted curves correspond to the BCS and BEC borders of 
the crossover region, respectively, as determined in Ref. 1421 . 
The potential strength V is in units of the critical potential 
strength V c = 47r/(mfco). 



(27r, 1/tt) of the parameter fc^pair were taken in Ref. Ha 
as representative of the BCS and BEC borders of the 
crossover region, essentially on the basis of the behaviour 
of the chemical potential as a function of the parameter 
fc^Cpair itself. 

We can see from Fig. 2 that the curve corresponding 
to the peak position V p extracted from the fidelity sus- 
ceptibility lies as expected in the middle of the crossover 
region for kp/ko < 1. For higher densities it approaches 
instead the BCS edge of this region. This may be due 
to the fact that at these densities and for fc^£ pa ; r of or- 
der unity (as in the middle of the crossover region), the 
range of the attractive potential becomes larger than the 
Cooper pair size. This favours the clustering of pairs 
(and eventually leads to an instability for sufficiently 
strong attraction, because of dominant pair-pair attrac- 
tive interactio n 24 i 43 i 44 ) . The pair correlation length will 
thus be increasingly influenced by inter-pair correlations 
rather than by intra-pair correlations. This implies that 
at high densities £ pa i r tend to overestimate the actual 
radius of a Cooper pair. At high densities the curve ob- 
tained from the peak position in the fidelity susceptibil- 
ity, which is not based on £ pa i r and is thus not influenced 
by this effect, is thus arguably a better indicator for the 
position of the crossover region than what is obtained 
from the calculation of fc_F^ pa ir itself. 

Note finally the merging of the three curves into the 
single point V = V c for vanishing kp/ko. This is ex- 
pected on physical grounds since in the two-body prob- 
lem a qualitative change in the ground state (from a de- 
localized to a localized wave function) occurs precisely 
at V = V c . It is indeed easy to verify that the fidelity 



FIG. 3: Fidelity susceptibility (normalized for convenience to 
the peak value) for the three-dimensional attractive Hubbard 
model at filling values n = 0.25, 0.05, 0.005 from left to right. 



susceptibility calculated over the two body wave function 
diverges at V — V c . At low densities, the fidelity suscep- 
tibility approaches the two-body behaviour and is thus 
peaked aroud V c . More generally, when the interparti- 
cle distance gets larger than the range of the potential 
the universal behaviour described by the contact poten- 
tial is recovered. Since, as we have seen above, for this 
potential the crossover occurs in the region \kpap | > I, 
it is then clear that when kp is vanishing the scattering 
length ap is bound to diverge in order to keep the prod- 
uct |fci?a_F| 1, thus pinning the crossover region close 
to V — V c where the scattering length diverges. 



IV. ATTRACTIVE HUBBARD MODEL 

We analyze finally the fidelity susceptibility calculated 
over the BCS wave function for the attractive Hubbard 
model with nearest-neigbour hopping: 



H = —t 



E 



U cJ t CiTc| 



(13) 



where indicates a sum over nearest-neighbor pairs 
and c| CT creates one electron with spin a in the Wannier 
state centered around the lattice site i. For an s-wave 
gap, which does not depend on k, the wave vector sums 
in the gap and particle number Eqs. (O,® and in Eq. © 
for the fidelity susceptibility are more efficiently calcu- 
lated by converting them into integrals over the energy 
and by using the density of states appropriate for the 
lattice kinetic energy dispersion. 

The resulting fidelity susceptibility for the three- 
dimensional attractive Hubbard model is shown in Fig. 3, 
at three different values of the filling factor n (parti- 
cle number per lattice site). For decreasing filling, the 
peak in the fidelity susceptibility progressively narrorws 
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FIG. 4: Attractive Hubbard model: position of the peak in 
the fidelity susceptibility x{U) in the plane (U/t,n) in three 
(a) and two (b) dimension. In two dimensions the curves 
corresponding to the BCS and BEC borders of the crossover 
region as determined in Ref. [4^ are also reported for compar- 
ison (dashed and dash-dotted curves, respectively). 



in width, while its position approaches the critical value 
(U c ~ 7.9t) above which a bound state appears in the 
two-body problem. The width of the crossover region as 
extracted from the half-height width in the fidelity sus- 
ceptibility thus shrinks for decreasing density. This over- 
all behaviour of the fidelity susceptibility suggests a pro- 
gressive approaching of the BCS-BEC crossover towards 
a quantum phase transition for decreasing density. It is 
in fact the presence of a finite density in the many-body 
problem which smears out over a finite coupling range 
the sudden change occurring at U c in the two-body prob- 
lem. In addition, dynamical mean field calculations^^ 
have shown the occurrence of a quantum phase transition 
between a Fermi liquid state and a paired state when su- 
perconductivity is artifically suppressed in the attractive 
Hubbard model. This transition, which also approaches 
the critical value U c for the two-body problem for van- 
ishing density, is transformed into a crossover when su- 
perconducting correlations are restored (in the very same 
way the metal-insulator transition is transformed into a 



crossover within the antiferromagnetic state in the half- 
filled repulsive Hubbard model). The behaviour of the 
fidelity susceptibility at low densities can thus be inter- 
preted as the remnant, within the superconducting state, 
of such underlying quantum phase transition in the nor- 
mal state. 

Fig. 4 compares finally the peak position in the fidelity 
susceptibility x{U) in three and two dimensions (full lines 
in panel (a) and (b), respectively). In two dimensions 
the curves corresponding to the conditions fcF^pah- = 27r 
(dashed line) and fcf^pair = 1/tt (dash-dotted line), ob- 
tained previously in Ref. [42| are also reported (the three 
dimensional Hubbard model was instead not considered 
in Ref. l42l ). From Fig. 4 (b) we can see that at low 
filling the peak position in x(U) nes m the middle of 
the crossover region, while at higher filling it bends to- 
wards the BCS border of the crossover region, similarly 
to what found for the three-dimensional NSR continuum 
potential. We think that also in this case the increasing 
importance of inter-pair correlations at higher densities 
explains the relative behaviour of the two curves in this 
regime. Note, in this respect, that for the Hubbard model 
it is the lattice spacing to provide the additional length 
scale relevant at high densities, playing the same role as 
the finite range of the interaction for the NSR potential. 

Note, finally, that in two dimensions the three curves 
reported in Fig. 4 (b) tend to the "critical" value for 
the two-body problem (U = 0) only at extremely low 
fillings (recall that in two dimensions, a bound state oc- 
curs in the two-body problem as soon as the attraction 
U is switched on), while in three-dimensions (Fig. 4 (a)) 
the peak position in the fidelity approaches the two-body 
critical value U c ~ 7.9t more progressively. The be- 
haviour of the three curves at low fillings can be explained 
by recalling that in two dimensions the binding energy of 
the two body problem eo oc exp(— t/U), while at low den- 
sities the Fermi energy relative to the bottom of the band 
is proportional to the filling (cf ccn), such that when eo 
and tF are of the same order (as it occurs for a crossover 
curve) one has n oc exp(— t/U). This explains the rapid 
collapsing to zero of the three curves when U/t < I . 



V. CONCLUDING REMARKS 

In this paper, we have tested the usefulness of the fi- 
delity susceptibility in the context of a crossover problem. 
We have considered specifically the BCS-BEC crossover, 
for which the BCS wave function is known to provide a re- 
liable description of the ground-state properties, and cal- 
culated the fidelity susceptibility over the BCS wave func- 
tion for several models exhibiting the above crossover. 
For the three-dimensional contact potential, we have 
found that the peak position in the fidelity susceptibil- 
ity and its width in terms of the dimensionless coupling 
parameter (/cfOf) -1 are in full agreement with previous 
knowledge about the position and range of the crossover 
region in this model (which has been widely studied in 
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the literature). 

For the finite-range NSR potential, the curve resulting 
from the peak position in the fidelity susceptibility was 
compared with the crossover "phase diagram" previously 
obtained in Ref. where the borders of the crossover 
region were defined on the basis of the value of the ratio 
(fcp^pair) between the pair correlation length and the av- 
erage interparticle distance. For densities such that the 
average interparticle distance does not get smaller than 
the range of the potential, the peak in the fidelity lies as 
expected in the middle of the crossover region, while for 
higher densities it tends to approach the BCS border of 
the crossover, as defined from the value of kp £ pa i r . We 
have attributed this behaviour to the increasing impor- 
tance of inter-pair correlations in the extraction of £ pa ir 
from the pair correlation function at high densities. 

We have considered finally the attractive Hubbard 
model. The results in two dimensions compared favor- 
ably with the corresponding crossover "phase diagram" 
of Ref. [42. In three dimensions, we have argued that 
the fidelity susceptibility is able to evidence within the 
superconducting state the traces of an underlying quan- 
tum phase transition, which has been found previously 
for the normal state when the superconducting correla- 
tions were artifically suppresse d 45 ' 46 . 

In summary, from our analysis we conclude that the fi- 
delity susceptibility provides a useful tool to characterize 



the width and position of the crossover region in crossover 
problems, which is especially appealing because of its 
"intrinsic" and quite general character. The definition 
and location of the crossover region based on the fidelity 
susceptibility do not depend, in fact, on specific quanti- 
ties like, for instance, the pair correlation length £ pa j r in 
the BCS-BEC crossover, which need to be changed when 
passing from one crossover problem to another one. Even 
though our analysis was based on the use of an approx- 
imate ground state wave function (the BCS trial wave 
function), the widely tested reliability of such a wave 
function in the context of the BCS-BEC crossover give 
us confidence in the robustness of our results. It will be 
however interesting to consider in future work alterna- 
tive methods, like e.g. Quantum Monte Carlo or Density 
Renormalization Group methods, or alternative crossover 
problems to confirm our main conclusions and place them 
in a broader context. 
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